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Abstract: The discovery potential of the LHC is investigated for the minimal 
anomaly-mediated supersymmetry breaking (mAMSB) scenario, using the ATLAS 
fast detector simulator, including track reconstruction and particle identification. 
Generic supersymmetry search cuts are used to map the 5 cr (and > 10 event) 
discovery contours in the m -m 3 /2 plane. With 100 fb _1 of integrated luminosity 
the search will reach up to 2.8 TeV in the squark mass and 2.1 TeV in the gluino 
mass. We generalise a kinematical variable and demonstrate that it is sensitive to the 
small chargino-LSP mass splitting characteristic of AMSB models. By identifying 
tracks from chargino decays we show that the Wino-like nature of the LSP can be 
determined for a wide range of chargino lifetimes. 
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1. Introduction 



Weak scale supersymmetry (SUSY) is the most promising solution to the technical 
hierarchy problem of the Standard Model (SM). However supersymmetric partners 
of known particles have not yet been discovered by experiment so if SUSY exists it 
must be broken. Dynamical supersymmetry breaking is usually assumed to occur 
in a 'hidden sector' which is isolated from ordinary particles and interactions. It 
is then communicated to the fields of the minimal supersymmetric standard model 
(MSSM) by some messenger (s). This must be achieved while maintaining both the 
theoretically appealing aspects such as naturalness [1,2] and respecting experimental 
constraints on masses, flavour-changing neutral currents and CP violation. 

SUSY breaking must be communicated to the observed particles by some inter- 
action felt by both the visible and hidden sectors. The most commonly discussed 
methods have been flavour-blind gravitational (SUGRA), or gauge (GMSB) interac- 
tions. 

More recently an alternative mechanism known as anomaly-mediation (AMSB) 
has been proposed [3,4] in which a conformal anomaly in the auxiliary field of the 
supergravity multiplet transmits SUSY-breaking to the observable sector. There will 
be an anomaly-mediated contribution to the gaugino masses in any hidden-sector 
model, but where no other direct contribution is present AMSB will be the leading 
effect. 

Anomaly mediation provides a potential solution to the SUSY flavour problem, 
in a highly predictive model. One undesirable feature of pure anomaly-mediation is 
that the slepton has a negative mass-squared. There are various ways to solve this 
problem [3,5-16], the simplest of which is the addition of a universal scalar mass at 
the GUT scale. The unmeasured parameters of the minimal model (mAMSB) are 
then m 3 / 2 - the gravitino mass; m - the universal scalar mass; tan/5 - the ratio of 
the vacuum expectation value of the Higgs fields; and the sign of the \i parameter 
multiplying HiH 2 in the superpotential, the magnitude of which is fixed from the 
condition of correct electroweak symmetry breaking. 

The sparticle spectra for mAMSB have been calculated in [17-20]. With pure 
anomaly-mediation, (i.e. m = 0) the gaugino masses are proportional to their beta 
functions: 

M, = ^im 3 / 2 (1.1) 

9i 

where gi are the gauge coupling constants with i = 1,2,3 indicating the gauge 
group, j3 gi are their corresponding renormalisation group beta-functions, and m 3 / 2 
is the gravitino mass. AMSB therefore predicts that the gaugino masses are in the 
approximate ratios M x : M 2 : M 3 3 : 1 : 7 so that the Wino (rather than the more 
conventional Bino) is the lightest supersymmetric particle (LSP), and the gluino is 
nearly an order of magnitude heavier than the LSP. 
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The prediction in AMSB of a Wino-like LSP has interesting phenomenological 
consequences. The most striking of these is that the lightest chargino is nearly mass- 
degenerate with the lightest neutralino. Near-degenerate particles are not unusual 
in SUSY phenomenology, but with AMSB one of these particles is the LSP. Since 
R-parity is assumed to be conserved in our analysis (unlike e.g. [21,22]) the xt ma Y 
only decay into the Xi so the small mass difference 

AM^=m(xf)-m(x?) (1.2) 

means that the lightest chargino may have a lifetime long enough to be detected at 
collider experiments. 

In the present paper we examine the AMSB discovery potential of the LHC 
general-purpose detectors from several points of view. We concentrate on the AT- 
LAS detector but similar considerations will also apply to CMS. In the next section 
we review the properties of the recently proposed benchmark points for minimal 
anomaly mediation. In sec. ^ we highlight some of the important features of the 
ATLAS detector. We explore the potential reach of generic SUSY search techniques 
in detecting minimal anomaly mediated supersymmetry in sec. |j. While the generic 
search reach is dependent on the mAMSB spectrum, the same is not true for the 
classic anomaly-mediated signature - the Wino-like LSP - which will apply beyond 
the minimal model. The signatures for identifying Wino-like LSPs at the LHC are 
explored in sec. Finally in sec. |6] we examine constraints on mAMSB from mea- 
surements of the cosmological relic density, the muon anomalous magnetic moment, 
and the branching ratio B — > X s ~/. 

2. Benchmark points 

A set of benchmark points and 'slopes' or model lines was suggested for study at 
Snowmass [23]. Of the eleven points, only one (SPS 9) applies to mAMSB. That 
point has the parameters: mo = 450 GeV, m 3 / 2 = 60 TeV, tan/5 = 10, with \i > 0, 
and lies on the model line "slope" itlq = 0.0075 x m 3 / 2 , where m 3 / 2 can vary. Part 
of the sparticle spectrum for SPS 9 is shown in fig. [I]. The lightest sparticles - 
the xt an d Xi — have masses about 170 GeV while squark and gluino masses are 
about 1.25 TeV, so one would expect copious sparticle production at the LHC. The 
Wino-like character of the lightest chargino and neutralino increases the relative 
cross-section to non-coloured sparticles as compared to e.g. SUGRA models. Indeed 
the HERWIG-6 . 3 [24-26] Monte-Carlo event generator gives the inclusive SUSY cross- 
section as 3.9 pb, of which about 0.5 pb is to squarks and gluinos. The chain 
Q X2 ~* I ~^ Xi is available, 1 and has a large branching ratio, so good information 
about the squark, slepton, and at least two neutralino masses could be extracted 
from kinematic edges in the 11, Iq and llq invariant masses [27-29]. 

1 Note that SM particles are omitted from decays when this can be done without ambiguity. 
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Figure 1: Part of the sparticle spectrum at the Snowmass point SPS 9 which has m,Q = 
450 GeV, m 3 / 2 = 60 TeV and tan/3 = 10, (j, > 0. Solid black lines indicate branching ratios 
(BR's) greater than 10%, dashed blue lines show BR's in the range 1% — > 10%, while red 
dotted lines show BR's in the range 0.1 — > 1%. The sparticles are displaced horizontally 
for clarity. 



The en and e~i masses are approximately equal (fig. |1|) because in mAMSB the 
slepton masses are principally determined by m . A more detailed analysis would 
be required to determine if the and masses could be separately measured from 
two nearby edges in the 11 invariant mass distribution, as was noted in [29] for a 
similar point. The decay of the right-sleptons to the Wino-like LSP is suppressed, 
but not sufficiently to produce a measurable displaced vertex from slepton decay 
(r ~ 10~ 16 s). In fact even if the Bino-like component of the LSP becomes extremely 
small, the decay e.R — > f\ — > x? wm remain unsupressed. This is particularly true at 
high tan (3 since the mixing between the left- and right-handed staus increases with 
the tau Yukawa coupling. 

Because the gluino has a similar mass to the heavier squarks at SPS 9 it decays 
primarily to ti and bi, meaning that gluino production will lead to large numbers 
of 6-quarks. This is not a general (or unique) feature of mAMSB, and while vertex 
tagging of jets could help distinguish this particular point from the SM it is not used 
in our analysis. 
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Figure 2: Part of the sparticle spectrum at the Point d'Aix: m,Q = 1000 GeV, m 3 / 2 = 
30 TeV, tan/3 = 30, \x > 0. Branching ratios as for fig. [I]. The sparticles are displaced 
horizontally for clarity. 



Another set of benchmarks, the 'Points d'Aix' [30] chooses the mAMSB point 
mo = 1000 GeV, m 3 / 2 = 30 TeV also with \i > 0, and with the higher value tan (3 =30. 
With a xt mass of 85 GeV, this point violates the limits from the LEP search [31] 
for e + e~ — > XiXil-i which places a lower limit on the xt mass of about 92 GeV for 
almost all values of the mass difference AM^ . 

The mAMSB Point d'Aix, despite being ruled out, still has a number of features 
which are of qualitative interest (fig. The light gluino means that the decays 
g — > ti and bb are no longer kinematically allowed. This means that the major 
source of b quarks is from X234 Xih followed by h — > bb, so fewer heavy quarks are 
produced than for the Snowmass point. The heavy sleptons can no longer participate 
in the various chains q — > Xx ~^ I ~^ X° y so slepton mass measurements would be 
extremely difficult at a hadron collider. There will, however, be large numbers of 
events in which xt s are produced, either directly or from g — > xtlQ^ so a signature 
based on the near-degeneracy of the xt an d Xi is robust. 
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3. The ATLAS detector 

ATLAS is one of two general-purpose experiments which will explore a wide variety 
of high pt physics with proton on proton collisions at y/s = 14 TeV the LHC. The 
LHC design luminosity is (10 34 cm _2 s _1 ), with a bunch-crossing period of 25 ns. 

The detector design, performance and physics potential are described in [27]. 
The main components are: 

• An inner tracking detector contained within a 2 Tesla solenoidal field which 
comprises three different elements, each of which has acceptance up to \rj\ < 2.5 
and 2ir in 0, where rj = — log[tan(0/2)] is the pseudorapidity, 9 is the polar 
angle from one of the beams, and is the azimuth. In order of increasing 
radius in the transverse direction (r in mm) the sub-detector components are: 
a pixel detector (43 < r < 132), a silicon microstrip detector (300 < r < 520) 
and a straw-tube transition radiation tracker with electron-identification ability 
(560 < r < 1070). 

• A highly segmented lead/liquid- argon electromagnetic calorimeter allows elec- 
tron and photon identification. The hadronic calorimeter employs steel tiles 
with a scintillator in the barrel region < 1.7) with granularity Arj x A0 = 
0.1 x 0.1, and copper/liquid-argon in the endcaps (1.5 < \r)\ < 3.2). A dense 
forward calorimeter with granularity 0.2 x 0.2 extends coverage to \r)\ < 5. 

• A muon system consisting of both precision and fast (trigger) detection el- 
ements. It is constructed within a toroidal magnet system and extends to 
transverse radii of about 10 m. 

4. LHC reach for mAMSB 

While some of the phenomenological features are particular to AMSB, one might 
expect to be able to distinguish AMSB from the SM using the same types of cuts - 
based on leptons, jets and missing energy - that are applied in mSUGRA analyses. 

The sensitivity of the LHC to mAMSB has been demonstrated [29] for one point 
(m = 200 GeV, m 3 / 2 = 35 TeV, tan/5 = 3, fi > 0) with relatively light sparticles, 
where the sparticle spectrum was investigated in detail. In [32-34], the signatures 
for AMSB at a future linear e + e~ or e~7 collider were investigated. In [35] the reach 
of the LHC was investigated using a simple generic detector simulation for 10 fb _1 
of integrated luminosity. The production of charged and neutral Winos via vector 
boson fusion was studied for AMSB in [36]. In that paper the LHC's reach was 
investigated for a signature consisting of two jets widely separated in pseudorapidity 
in association with missing transverse momentum. 
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In this section our the aim is to determine the reach of the LHC with a realis- 
tic detector simulator, using optimised but generic SUSY cuts and for 100 ftT 1 of 
integrated luminosity. 

4.1 Event simulation 

The mAMSB spectra were generated using ISAJET-7 . 63 [40,41] on a grid 100 GeV x 
5 TeV in the (mo, 7713/2) plane. In all cases the ratio of the Higgs vacuum expectation 
values (tan (3) was set equal to 10 and the sign of /x was positive. The mass of the 
top quark (important for electroweak symmetry breaking) was taken to be 175 GeV 
throughout. The dependence of some of the key sparticle masses on the input values 
of mo and 7713/2 is shown in fig. |3|. 

In the xf — > Xi decays, ISAJET does not include masses for the leptons, and 
does not contain multi-pion decay modes. Since the mass difference AM^ can be 
of the order of the mass of the muon, the lepton mass effects can be important in 
AMSB. To improve accuracy, the chargino decay modes calculated in [37,38] were 
implemented with pion form factors from [39] and massive leptons. The resulting 
chargino lifetime and branching ratios are shown in fig. |j. 

HERWIG-6 . 3 was used to produce 200 fb _1 of unweighted inclusive supersym- 
metry events for each point (with a minimum of 5 x 10 4 , and up to a maximum of 
5 x 10 5 events). 2 This was then scaled to give the expected number of events for J C = 
100 fb _1 , which corresponds to the first year of 'high luminosity' (10 34 cm _2 s _1 ) run- 
ning of the LHC. Background samples were generated with HERWIG for the production 
of W^+jets, Z°+jets, tt, and QCD 2^2 (excluding ti). For the W^+jets sample, 
the cross-sections were multiplied by the factor 

2(iV jet -l) 

(4.1) 

where A/j et is the number of jets, and p^ hr is the jet transverse energy threshold 
which was set to 10 GeV. This correction brings the HERWIG cross-section into better 
agreement with tree- level matrix element calculations [44] . A total of over 20 million 
background events were generated in logarithmic intervals in the HERWIG parameters 
PTMIN and PTMAX from to 7000 GeV. The background cross-sections are shown 
before and after preselection cuts in fig. || 

The events were passed through the ATLAS fast detector simulator, ATLFAST-2 . 50 
[45], which gives a parameterised detector response based on GEANT3 Monte-Carlo 
simulations [46] . Jets were found using the ATLFAST cone algorithm with a cone size 
AR = 0.4, and a minimum of 10 GeV. The loss in resolution from pile-up was 
simulated. Calorimeter cells with Et deposits below 1 GeV were not included in the 

2 No significant differences were found on comparing with HERWIG-6. 4 [42] which includes spin 
correlations as described in [43]. 
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Figure 3: The contours show the (a) gluino, (b) up left squark, and (c) lightest neutralino 
mass in GeV, and (d) the Am = xt ~ Xi ma ss difference in MeV as a function of mo 
and m-3/2- The other parameters are: tan/3 = 10 and fx > 0. The Snowmass point 
(mo = 450 GeV, = 60 TeV) is marked as a pair of concentric circles, through which 
passes the Snowmass slope "model line G" (dashed line) . The solid red regions are excluded 
because of lack of electroweak symmetry breaking (bottom right) or because of charge- or 
colour-breaking minima, or non-x? LSP (left). Am is greatest where [i is small, which 
occurs near to the region where electroweak symmetry is unbroken. 
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(a) (b) 

Figure 4: The (a) lifetime and (b) branching ratios of the lightest chargino as a function 
of the mass difference AM X = M(xt) ~ M(xi)- The ra pid decrease in the lifetime occurs 
at AM X = m w + where the single pion mode becomes available. The discontinuity at 
AM X = 1.4 GeV comes from the switch in the calculation from hadronic to partonic 
decay widths. The leptonic channels implicitly include the corresponding neutrino. The 
branching ratio to 7r + 7T°7r° was assumed to be equal to 7r + 7r~7r + . After [37-39]. 
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IpT calculation in order to more accurately model the expected resolution. The cell 
sizes were At] x A0 = 0.1x0.1 for \r)\ < 3 and 0.2x0.2 for 3 < \r]\ < 5. Otherwise the 
default ATLFAST parameters were applied. The parameterised tracking simulation 
used is described in appendix |B|. 

4.2 Optimisation of cuts 

Cuts were applied to leptons (electrons and muons only), jets, and missing transverse 
energy in a similar manner to those used in [27,35,47-49]. The variables to which 
cuts were applied were: 

(1) JfjTi missing transverse energy; 

(2) Pt(ji), Pt of the hardest jet; 

(3) Pt(J2), Pt of the next-to-hardest jet; 

(4) YIiPti scalar sum of the pt of jets in the event; 

(5) A/jet, number of jets in the event; 

(6) St, transverse sphericity (circularity) of the event; 

(7) A0(j!), difference in azimuth between hardest jet and $t vector; 

(8) Pt(I\)i Pt of hardest lepton (if any); 

(9) A0(£l), difference in azimuth between hardest lepton (if any) and jzJt vector; 

(10) Mt = a/2 PT(ei) $t (1 — cos(A0^!))), transverse mass of lepton and missing 
energy. Applied to single-lepton channel only, to reduce SM leptonic back- 
ground. 

The preselection cuts which were applied to cut down the background, par- 
ticularly at low pt, are shown in the third column of table [l]. The reduction in the 
background cross-section can be seen in fig. [5]. Hard preselection cuts were not placed 
on isolated electrons or muons since it is foreseen that one jet with pr > 100 GeV 
together with I^t > 100 GeV will be sufficient to provide the on-line trigger. Other 
trigger strategies which would improve the reach include: 

• Allowing a hard (~ 20 GeV) isolated lepton as a trigger could improve the 
selection for points for relatively low m in which the cascade decay q — > x° ~~ * 

^ s av& il a ble; 

• Track vertexing which in ATLAS can be applied at the second trigger level 
could improve the selection for events containing heavy quarks; 
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Table 1: Allowed values for each of the cuts described in the text. The preselection cuts 
(highlighted in bold italic) are shown in the third column (applies to cuts 1 to 5 only). 
The units of the variables (1-4,8,10) are GeV. 



• If all of the coloured sparticles are heavy then the dominant production of SUSY 
particles will be to the lightest gauginos. At leading order in perturbation 
theory, the production of x± s and/or xfs would be difficult to trigger on, but 
gaugino production in association with a high p? jet or photon and IpT might 
be observable [20,50]. 

A fraction of the simulated gaugino production events contain high jets from 
initial state parton showers. However the parton shower algorithm is based on the 
soft and co-linear approximations so this fraction will be smaller than that which 
would result from higher order matrix element calculations. Our simulation of the 
trigger can therefore be considered conservative. 

Different values of the cuts were applied to events with zero, one, two or three 
leptons (electrons and muons only) as reconstructed by ATLFAST, as well as to an 
inclusive lepton ('ptmiss') analysis. In the case of two-lepton events, different cuts 
were applied according to whether the leptons were of the same sign (SS) or opposite 
sign (OS). 

For each of these analyses the cuts were optimised in the order listed in table [L[ 
and allowed to take one of ten values shown. The significance of the signal was given 
by S/ \f~B where S and B are the number of signal and background events expected 
respectively for 100 fb _1 . For each variable, the cut was chosen to maximise the 
significance subject to the constraint that S > 10. The analysis was determined to 
be successful if significance greater than five was achieved with at least ten events 
passing for any set of cuts. 
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This technique will not necessarily generate the global maximum in S/ \J~B since 
the cut variables are not totally independent, but it is sufficient for the purposes of 
a large parameter-space scan, and it decreases the chances of over-fitting to a sparse 
background. 

4.3 Results 

The 5er (and > 10 event) discovery reaches for 100 ftT 1 of integrated luminosity in 
fig. ^ and fig. [F] show that the LHC will be able to distinguish mAMSB from the SM 
over a large range of parameter space. As expected, the results are similar to previous 
generic R-parity conserving SUSY searches [27,35,47-49] with reach extending up to 
squark masses of about 2800 GeV, or to gluino masses of about 2100 GeV, whichever 
is the lower. The discovery reach is similar to, but greater than that found in [36] 
for Wino production via vector boson fusion. 

In general the best reach is obtained in the single-lepton and the inclusive chan- 
nel. The same-sign dilepton signal is competitive at low m , when < m^o - 
a region which includes the model line SPS G - but this sensitivity decreases as 
m and m| increase. Points with large mo have a spectrum qualitatively similar to 
that in fig. ^|, with heavy squarks and sleptons. Such points can still produce lep- 
tons, mostly from heavy chargino or neutralino decays. Another source of leptons 
is from gluino decays to higgsinos along with a quark and anti-quark from the third 
generation, followed by leptonic t or b decay. 

Along the Snowmass model line 'G' (mo = 0.0075 x 777,3/2) the LHC will be able to 
measure discrepancies from the Standard Model up to m =1050 GeV, m 3/ / 2 =140 TeV, 
at which point the gluino mass is 2.76 TeV and m^ L = 2.70 TeV. The Snowmass 
point SPS 9 lies well within the discovery region for all of the six different analyses. 
The optimised cuts and the resultant significance are shown in table and |3|. The 
high significance implies that it will be possible to extract further information from 
the data, for example to determine sparticle masses. 

5. Distinguishing Wino-like LSPs 

Most of the mass-sum rules [19] for mAMSB involve terms like Mj — M\ L which 
would be extremely difficult to determine experimentally. The exception is the acci- 
dental near-degeneracy of the and e#, which might be observable at some points, 
as was noted in sec. However points with large m have relatively heavy sleptons 
which would have a small production cross-section at a hadron collider. The Wino- 
like LSP signature of AMSB, and the resultant near-degeneracy of the \t an d Xi 
will be the robust 'smoking gun' for anomaly mediation, and is applicable beyond 
the minimal model. 

The LHC, like LEP [31], or the Tevatron [20,50] or a future linear e + e~ collider 
[51] can use a variety of analyses in the search for long-lived charginos, each tailored 
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Figure 5: The cross-sections for the tt, + Jets, and Z® + Jets and QCD (excluding tt) 
backgrounds plotted as a function of the transverse momentum variable which is compared 
in HERWIG to PTMIN and PTMAX. The production cross-sections are denoted with dashed 
lines, while the cross-section to pass the preselection cuts are solid lines. Where no events 
passed the selection cuts (at low PT) arrows indicate 90% Poisson confidence limits on the 
selected cross-sections. 
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Table 2: The value of the cuts which maximised the significance (S/\/~B) for each of the 
analyses at the Snowmass mAMSB Point SPS 9 for fC = 100 fb -1 . The M T cut applied 
only to the single lepton channel and was kept at 100 GeV. 
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Figure 6: The 5 a (and > 10 event) discovery region for mAMSB in the mQ~m^/2 plane, 
is shown by the light blue diamonds, for the (a) no lepton, (b) single lepton, (c) two 
opposite-sign leptons, (d) two same-sign leptons channels. The other parameters are 
tan/3 = 10 and [i > in all cases. The solid red regions are excluded because of lack 
of electroweak symmetry breaking (bottom right) or because of charge- or colour-breaking 
minima or non-x? LSP (left). The low m 3 / 2 region excluded by the LEP limit [31] on the 
chargino mass is indicated by purple crosses. The Snowmass mAMSB point SPS 9 at 
tuq = 450 GeV, m 3/ / 2 = 60 TeV is marked with a pair of concentric circles, through which 
passes the Snowmass slope 'model line G' (dashed line). Contours of ul and g iso-mass 
are shown as solid lines. 
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M / GeV M / GeV 



(a) (b) 

Figure 7: As fig. || but for the (a) the three lepton channel and (b) the inclusive ('missing 
pt') channel. 
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B 
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S 


B 


s/Vb 


ptmiss 


3300 


151 


270 


2650 


0.25 


5200 


Lep 


950 


21 


206 


7800 


4.46 


3700 


1 Lep 


151 


0.00005 


21000 


930 


0.00005 


10 5 


2 Lep OS 


15 


0.00005 


2100 


62 


0.00005 


8800 


2 Lep SS 


13 


0.0005 


1900 


67 


0.0048 


970 


3 Lep 


206 


0.0001 


20400 


19 


0.00010 


1900 



Table 3: The expected number of signal (S), and background {B) events and the signifi- 
cance (S/y/B) for for each of the analyses, at the Snowmass mAMSB Point SPS 9 and the 
Point d'Aix, for J C = 100 fb _1 . For all six analyses and for both points the backgrounds 
are dominated by high vector boson production in association with jets. 
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M / GeV 

Figure 8: The number of X\ ~* X? decays expected within the central region of the 
detector (\rj\ < 2) with transverse decay vertices between 100 mm and 800 mm from the 
interaction point for integrated luminosity of 100 fb~ . The initial track from the chargino 
is required have pr > 10 GeV. All events are required to have passed the preselection cuts 
(table [l]) which ensure that the event will be triggered. The 10 event and 20 event contours 
have been smoothed for clarity. AAfcj is less than 200 MeV above and to the left of the 
dashed blue line on the bottom right-hand side of the plot. 



16 



to a different regime in AM^ , according to the xt lifetime. The various regimes are 
ordered by decreasing xt lifetime below: 

(1) AM^, < 771^+. If the pion decay mode is not available, the xt will ^e long lived 
(cr >lm) and can leave a track through the muon chambers. Analysis of this 
type of signature was performed in [52] for a GMSB model with long-lived f\ 
NLSP, and in [53] for an intermediate scale model with heavy stable leptons. 
The mass of the chargino can be measured using the muon detector as a time- 
of-flight system. Additional discrimination can be obtained by considering ^ 
information from the transition radiation tracker. Higher-mass sparticles can 
then be reconstructed by forming the invariant mass of the stable particle with 
jets and leptons. Provided the xt lifetime, although long, is short enough that 
a reasonable number of xt s decay within the inner detector, the mass differ- 
ence AM^j could also be determined from that fraction using the techniques 
described below. 

(2) 171^+ < AM^, < 200 MeV. This is the regime in which high px chargino tracks 
often decay within the body of the inner tracker to soft pions or leptons along 
with large amounts of missing transverse energy. The details of the detector 
resolution and track reconstruction algorithm are beyond the scope of this 
study, but such tracks should provide a striking signature provided they occur 
in sufficient numbers. The number of xt s which would be produced at the 
LHC and decay within a fiducial volume in the active material of the ATLAS 
tracker is plotted in fig. |as a function of itlq and m 3 / 2 . 

One-prong chargino decays such as xt ~ * Xi 7l+ an d xt ~^ Xi e+h, e will produce 
'kinks' at the point where the charged SM particle is softly emitted. In some 
cases the SM particle will go undetected and the high p? chargino track will 
seem to disappear. Such decays, producing track stubs, could also be detected 
with a dedicated off-line analysis. 

Even though rather strict triggering requirements have been made, there are at 
least ten such events for 100 fb" 1 over almost the entire plotted parameter-space 
provided m 3 / 2 ^ 185 TeV. This signature should be essentially free of physics 
background, so 'detector' backgrounds will dominate. Since ATLAS should 
have three barrel pixel layers with r < 122 mm with noise occupancy of less 
than 10 -5 , and average physics occupancy of the order of 10~ 4 the instrumental 
backgrounds should also be under control [54,55]. Such a search could therefore 
achieve a large reach and would be sufficient to identify a near-degenerate LSP 
model. 

(3) 200 MeV < < a few GeV. Chargino decays which occur before the 
first tracking layer can still be useful provided that the soft track from the SM 
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particle can be found, using e.g. electron identification, track impact parameter, 
isolation from other tracks, and the direction of the jziy 2-vector. 



The maximum value of AM^ in minimal AMSB with tan/3 = 10 is less than 
about 200 MeV except close to the region of no electroweak symmetry breaking. This 
leads to a decay distance cr ^ a few centimeters, and corresponds to category (2) 
above. However, other models such as 0-11 string models [37,38,56] have Wino-like xf 
and Xi with AM^ typically of the order of a GeV. Such models fall into category (3), 
which is potentially the most difficult at a hadron collider because the short-lived 
X\S will not be directly observable, and their softly emitted SM daughters suffer from 



a large background from other low-momentum tracks. In sec. |5.1| we define a range 
of model points which allows us to study this regime of shorter lifetimes. Some of 



the methods which can enable the LHC to probe this region are explored in sec. [572 



5.1 Definition of points 

The tree-level neutralino mass matrix is: 
/ 



Mi — mzcpsw rnzspsw \ 

M 2 m z Cf3C W -mzspcw 

-m z cps w mzCpCw -// 



V 



(5.1) 



and the tree-level chargino mass matrix is: 



M 2 \i2m w sp\ 



(5.2) 



where sw = sin^ ; c w = cos 6^, sp = sin/3, cp = cos/3 and we use the conventions 
of [57]. The mass difference m~+ — m~+ = AM^ is highly suppressed at tree-level, 
so the leading 1-loop correction can be important. It takes the form [58]: 



Ai\% (1 - loop) 



a 2 M 2 
Air 



m w 
~M~ 2 



m z 
M 2 



+ 5s 



w 



(5.3) 



where F(a) = dx{2 + 2x) log[x 2 + (1 — x)a 2 ], and is included in ISAJET. 

We are interested in the collider phenomenology of long-lived charginos with 
various lifetimes. From (|5.1|) and (|5.2j ), the leading tree- level mass difference term 



AM-^ oc l//i 2 . In order to explore different regimes in AM^, we define new points 
based on SPS 9 and another minimal anomaly-mediated point with mo = 500 GeV, 
m 3/2 = 36 TeV, fi > and tan/3 = 10. We then decreased the value of the /i 
parameter at the electroweak scale to the values shown in table 01. 
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Point 




m xt 
(GeV) 


AM* 
(MeV) 


cr 
(mm) 


Al ^ Al e u e 


Al ^ A1P U V 


SPS 9 


864 


171 


164 


56 


2.0 % 


0.2 % 


SPS-300 


300 


165 


886 


0.11 


17.0 % 


15.9 % 


SPS-250 


250 


159 


1798 


0.004 


21.9 % 


21.5 % 


A 


533 


107 


181 


34 


2.0 % 


0.3 % 


A-250 


250 


101 


766 


0.20 


15.4 % 


13.9 % 


A-200 


200 


97 


1603 


0.007 


22.5 % 


22.2 % 



Table 4: The mass of the lightest chargino, its mass difference with the LSP, decay length 
and leptonic branching ratios for six points. SPS 9 is the Snowmass mAMSB point, from 
which we define SPS 300, and SPS 250 by adjusting fi at the electroweak scale to produce 
different AM^ . Point A is the mAMSB point with m = 500 GeV, m 3/2 = 36 TeV, fi > 
and tan/3 = 10, from which \i is adjusted at the electroweak scale to produce A-250 and 
A-200. 

In order to maintain consistent electroweak symmetry breaking when decreasing 
/i, the Higgs soft mass parameters, m 2 Hi and m 2 H2 , must allowed to vary, as can be 
seen from tree-level equation, 

2 raji ~m 2 H2 tan 2 (3 1 2 
P tan 2 (3-1 2 z 

The points with adjusted fi are then no longer in the minimal anomaly-mediated sce- 
nario, since m? H and m 2 H2 are not those which would be predicted from the mAMSB 
parameters mo and 777,3/2- Decreasing /i has the side-effect of decreasing the masses 
of the higgsinos, and can change the phenomenology somewhat - for example by 
opening chains such as q — > xtl ~^ XiQ^ . When followed by leptonic Z° decay 
these chains could, in principle, be used to further constrain the dynamics. However 
the only requirement we make for our analysis is that the decay q — > XiQ occurs, 
producing highly boosted xf pairs in association with jets. 

5.2 Identifying chargino decay products 

In this section we demonstrate how the low-pr tracks from the charged SM daugh- 
ters of xt decays can be identified and used to constrain AM^. To do so we must 
distinguish these tracks from the large number of other low-p^ tracks from the par- 
ticles generated in the proton remnant interactions, and from any pile-up events. To 
avoid having to deal with the contamination from multiple pile-up, in this section 
we simulate the initial three years of 'low luminosity' (10 34 cm~ 2 s _1 ) running of the 
LHC, which is expected to provide 30 ftT 1 of integrated luminosity with an average 
of one inelastic collision per bunch-crossing. The loss of resolution from pile-up is 
therefore not simulated in this section. 
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Figure 9: The transverse momentum distribution of reconstructed charged tracks for 
the point SPS-250 (defined in table g) after selecting events with $ T > 200 GeV. The 
thin black line shows all tracks; the distribution for particles originated from X\ decays 
(the signal) is highlighted by the thick red line, where the thickness shows the statistical 
uncertainty; those from the minimum bias event and from the soft underlying event are 
denoted by circles and asterisks respectively. 




Figure 10: The number of tracks surviving subsequent cuts, after a preselection cut 
$ T > 200 GeV for the point SPS-250 (defined in table g). The signal tracks (from xf 
decays) are also plotted separately, as are tracks from the HERWIG soft underlying event and 
the additional minimum bias QCD event. Not all cuts are used in any particular analysis. 
For this point the most effective cuts are those on the px of the track, and on selecting 
leptons ('lept'), while further cuts on the transverse impact parameter (do) and isolation 
(rj — <p) improve the selection. 
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The HERWIG soft underlying event model, based on the UA5 collaboration [59] pp 
Monte-Carlo, is known to underestimate the number of tracks with pt > 1 GeV [60], 
because it does not model the semi-hard physics of multiple parton interactions. To 
ensure that our analysis was insensitive to contamination from such tracks, a sample 
of 5000 inclusive QCD events with PTMIN = 4.5 GeV was generated. We refer to these 
as 'minimum bias events' because they simulate (better than the default HERWIG soft 
minimum bias option) the expected minimum bias distribution of tracks, including a 
tail at higher p? due to the onset of hard scattering. The particles from this minimum 
bias event were added to the SUSY event in addition to HERWIG's usual underlying 
event. 

An initial selection requiring a high pt jet and large missing transverse energy 
was made to trigger the event, reduce the SM background, and select those SUSY 
events in which the neutralino was highly boosted. In this section we are dealing 
with the short-lifetime regime (cr < 1 mm), so this boost will almost certainly not 
be sufficient to allow the chargino to live long enough to decay in the body of the 
detector. However a boosted chargino means that even particles which are very softly 
emitted in the chargino rest frame will have a transverse momentum 

Ptx ~ x m x 

Al 

as measured in the detector, where x is the charged SM daughter particle and p T ^o 
is of the order of the missing transverse momentum. This extra px assists the re- 
construction of the SM daughters of charginos and helps distinguish them from even 
softer particles coming from the underlying event (fig. ^). 

The tracks from SM daughters of chargino decays have other properties which can 
help in their identification. For larger AM^ the xt branching ratio to each of muon 
or electron and associated neutrino can be up to nearly twenty percent (fig. f|b). In 
the ATLAS experiment electrons and muons will be distinguished from other tracks 
by transistion radiation and by penetration respectively (see appendix |B|) . 

Much of the background is associated with heavy-quark decay and can be re- 
moved by applying isolation cuts (table ||) which require that for any track, i, to be 
a candidate: 



no other track is found with Ai?y = y/ (Ar/^) 2 + (A0jj) 2 less than some value 
AW 



jmin. 



the invariant mass of the track with another track, j, is greater than M™ m for 
all j, where it is assumed that m; = rrij = m n +; 

no other track with R±j < R imp has Oq j 'a(do ) > 3, where is the transverse 
impact parameter, and ct(g^) is its uncertainty; 



Si do / a (do ) < °"imp where the sum is over all j ^ i with Rij < R 



imp' 
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JT UlllLo 




cpcocf) A-9f)f) 


Event: 




500 


500 




min 

Pt(ji) 


400 


400 




Omin 


0.05 


0.05 


Track: 


,v-,min 
FT 


2.0 


2.0 




max 
FT 





1U 




Jmin 


0.03 







-rJHlllX 


0.3 


0.1 




1 ^ 1 max 


1.5 


2 




AR™ a 


0.45 


0.2 






0.9 


0.4 




Rimp 


0.4 


0.4 




C^imp 


3.0 


3.0 




Particle 


any 


£ E e, [i 



Table 5: The cuts applied to events and to tracks for the different points in table ||. The 
missing transverse energy ($r)> leading jet transverse momentum (pt(ji)), and transverse 
sphericity (St) cuts were applied to the whole event, while the other cuts were applied track 
by track. The track transverse momentum was required to lie in the range p™ in — * p™ ax ; 
its transverse impact parameter in the range ci™ 111 — ► ci™^; and the absolute value of its 
pseudorapidity was required to be less than |r/| max . The other cuts are described in the 
text. Energy, momentum and mass units are GeV; the impact parameter, do, is measured 
in mm. 
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m , v -m z / GeV m Tx -m x / GeV 

(a) (b) 

Figure 11: The tuta — m y0 distribution for the points (a) SPS-250, and (b) A-200 
(defined in table f|). These two points have similar values of AM^, and the same cuts 
were applied to both. Only events in which exactly two leptons passed the selection cuts 
are plotted. The signal tracks are those from xt an d X\ decays. Events where only one of 
the two tracks has been successfully identified are plotted separately. The vertical dotted 
line indicates AMy, = m~+ — m^o. 
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Figure 12: The my2 — distribution for (a) the point SPS-300, and (b) the point 
A-250. These two points have similar values of AMjj, and the same cuts were applied to 
both. Only events in which exactly two particles (of any type) passed the selection cuts 
are plotted. The peak is closer than in fig. [ll] to the upper kinematic limit at niT2 — m x<( = 
AM^ (dotted line) since there are only two missing particles - the two neutralinos. 
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5.3 Results 



We selected those events in which precisely two tracks satisfied the track cuts listed 
in table |5|. The numbers of tracks from chargino decays and from other sources 
which passed subsequent cuts are shown in fig. [LCj for the point SPS-250. We then 
considered the Cambridge rriTx variable which (as described in appendix remains 
sensitive to AM^ even when there are uncertainties in the neutralino mass or the 
missing transverse momentum, providing an important handle on Wino-LSP physics 
at hadron colliders. We assume that both the squark mass scale, rriq, and the Xi mass 
have been previously determined from the measurement of other kinematic edges as 
described in [27-29,49], and plot distributions of m TX — m^o which has the property: 

<(m TX -m 2 o) < (m-+-m^o) = AA%. (5.4) 
For the two points with AM^ « 1.7 GeV, we selected isolated leptons and 



plotted distributions of m*r 4 — m^o. As can be seen in fig. [n], A-200 has more signal 
events passing the cuts than SPS-250, but the peak position is the same in both 
cases. 

As AM^j decreases, the leptonic branching ratios decrease (see fig. ^b) and the 
detector's ability to distinguish the lower leptons from hadrons is diminished (see 
appendix ||). For this reason, at the two points with AM^ « 800 MeV we select 
events with two tracks of any type satisfying the cuts in the third column of table [5[ 
Most of these tracks will be pions because of the large branching ratio to 7r =b or tt^tt . 
The m T2 — m^o distributions for these two points are plotted in fig. [L2 . 



For all four points the narrow peaks indicate that AM^ is of the order of 1 GeV. 
This confirms the Wino-like nature of the LSP, providing the 'smoking gun' signature 
for anomaly mediation. 

Each rriTx distribution depends principally on AM^, on the neutralino mass 
and the momentum distribution of the charginos. The latter depends largely on rriq 
so if both niq and the lightest neutralino mass were already measured, then AM^ 
can be measured by fitting to each m TX — m^o distribution. 

To demonstrate that m*T4 can indeed make a quantitative measurement of AM^ , 
a phase-space Monte-Carlo program was used to generate very simple 'events' in 
which pairs of squarks decayed via the chain q — > xi 1 ~ * Xi e v e Q- m Ti — m x ( { 
distributions were produced for those events in which fa > 500 GeV, and the peak 
determined from a gaussian fit. The correlation between the fitted peak positions and 



the input values of AM^ is shown in fig. [L3|. The peak of the rriT4 — m^o distribution 



for point A-200 (fig. IT|b), was likewise determined from a gaussian fit. As can be 
seen from fig. [13| the mass difference AM^ can be measured at that point with a 
statistical uncertainty of approximately 150 MeV (10%). 
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Figure 13: The diagonal line shows how the tjita — m^o peak value depends on AM^. 1 . 
The thickness of the line indicates the uncertainty in the peak position from a gaussian fit 
to the distribution combined in quadrature with a 10% uncertainty in the LSP mass and 
a 10 % uncertainty in the squark mass scale. The peak of the tuta — m x° distribution for 
the point A-200 is marked with a star at its input value of AM^ ; the error bar shows its 
uncertainty. 



6. Other constraints 

6.1 Cosmological relic density 

If R-Parity is conserved then the lightest supersymmetric particle can be a good can- 
didate for the cold (non-relativistic) dark matter hypothesised by astrophysicists and 
cosmologists. The cold dark matter contribution to critical density of the universe 
provided by the LSP is given by: 

<^CDM = ^CDM^- = 1 ~7T^ 

Pcrit/h 2 

where p cr i t = 3H 2 /8ttG is the critical density, n^o is the number density of LSPs, 
and H = h x 100 km s" 1 Mpc" 1 is the Hubble constant. 

Astronomical estimates from measurements of the acoustic power spectrum of the 
cosmic microwave background anisotropy [61-65] suggest that ujqdm = 0.106±0.010. 
Further constraints from the 2dF Galactic Redshift Survey [66] lead to c^cdm 
0.1 1 I ± 0.0091, assuming a flat universe (Ejf2; = 1). 
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In AMSB the Wino-like xi an d Xi undergo rapid annihilation though reactions 
such as Xi Xi ~ * W + W~ and Xi Xi ~ > Q q'- The relic density of cold dark matter 
was calculated for AMSB using the program micrOmegas [67], which includes the 
annihilation channels listed in table ^|. The results, in fig. |T4], show that the value 
of cg>cdm is small throughout the region, and about a factor of two less than was 
estimated in [4]: 



no-coanh _ g x 1Q -4 x f X 1 



cdm -VlOOGeV, 
which neglected the co-annihilation channels. 

Thus in AMSB the LSP does not give any problems with relic overabundance, 
but suffers the opposite problem - that sparticles in thermal equilibrium in the early 
universe will not produce sufficient neutralino dark matter. It has been suggested 
in [68] that decays of cosmological moduli - which are low mass particles predicted 
by string theory when supersymmetry is broken - could have produced Winos with 
sufficient abundance to be of astrophysical interest. 

6.2 Muon g-2 

It is well know that the gyromagnetic moment of the muon, can sensitive to 
sparticle interactions through loop corrections [69-73]. The recent BNL measure- 
ments [74,75] dominate the world experimental average of the positive muon anoma- 
lous moment, 

a™ L = - 2)/2 = 11 659 203(8) x 1(T 10 . 

Standard Model calculations of a M have been reviewed in [76,77]. Recent calculations 
[78,79], which include a correction to the sign of the pion pole part of the hadronic 
light-by-light contribution [80-84], were combined in [75] to give: 

11 659 177(7) x 10~ 10 < af 1 < 11 659 186(8) x 10~ 10 . 

The SUSY contribution is typically dominated by loops involving charginos and 
neutralinos [76,85,86], and is proportional to tan f3 in the high tan (3 limit. The AMSB 
contributions to have been calculated in [18,86]. In fig. |TJ we plot the a^ MSB for 
tan/? = 10 and 30, where again we have made use of the program micrOmegas. 3 
The AMSB contribution is small compared to both the experimental and theoretical 
uncertainties, except at small m , 7713/2- 

6.3 B -> A s7 

The inclusive radiative decay B — > X s 7 is sensitive to sparticle properties through 
radiative corrections involving charged higgs, chargino and t loops (see for example 
[18,87-89]). 

3 Note that ISASUSY produces AMSB points with negative Mi and Mi so the sign correlation 
between /j, and a^ USY is opposite to that in [86]. 



26 



Experimental measurements of the branching ratio from the CLEO [90], BELLE 
[91] and BaBar [92] collaborations give: 



B(B -> X s7 ) CLEO = 3.21 ± 0.43 ± 0.27+°;}^ x 10" 4 
B(B -> X s7 ) BELLE = 3.39 ± 0.53 ± 0.421°;^ x 10~ 4 , 
B(B -> X s7 ) BaBar = 3.88 ± 0.36 ± 0.37±g;^ x 10~ 4 , 

(6.1) 

where the errors are statistical, systematic, and from theory respectively. Renormal- 
isation scale uncertainties [93] lead to a 10% (±0.3 x 10 -4 ) uncertainty in the SM 
prediction. The mAMSB contributions were calculated using mi cr Omegas and are 



plotted in fig. [16|. The combined experimental and theoretical uncertainty does not 
allow us to constrain the mAMSB parameter space at better than 2a for either of 
tan/? = 10,30. 

The overall effect of the B —>■ X s7 and to a greater extent the g M — 2 constraint 
is to disfavour the low mo, m.3/2 region especially when tan/? is large. 



7. Conclusions 

If anomaly mediated supersymmetry is present at the 1 to 2 TeV scale, the LHC 
will observe excesses in various multi-lepton + fa channels. We have used generic 
supersymmetry search procedures, and a realistic detector simulation, to investigate 
the ability of the experiments at the LHC to discover AMSB scenarios. By selecting 
events with very large missing energy and identifying tracks from chargino decays, 
the Wino-like nature of the LSP can be determined, and m,-+ —m^o can be measured 
over a large range in parameter space. Careful study of the tracks left by \\ s an d 
their decay products can give clear evidence for this class of models even in cases 
where the chargino is shorter-lived than predicted in minimal AMSB. 
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Figure 14: The relic density of the lightest neutralino, f^cDM^ 2 with fi > 0; (a) for 
tan/3 = 10 and (b) for tan/3 = 30. The meaning of the symbols is explained in fig. |6|. 
The astrophysical limits on f^cDM^- 2 are discussed in the text. The contours are in units 
of 10" 3 . 
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XiXi - 
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15 % 


xixi - 


-> boson pair 


12 % 



Table 6: The main cosmological relic annihilation channels for the Snowmass point SPS 9. 
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Figure 15: The SUSY contribution to the anomalous magnetic moment of the muon, (in 
units of 10 -10 ) with \x > 0; (a) for tan/3 = 10 and (b) for tan/3 = 30. The meaning of the 
symbols is explained in fig. @. 
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Figure 16: The branching ratio B{B — ► X s j) with /j, > 0; (a) for tan (3 = 10 and (b) for 
tan /3 = 30. The contours are in units of 10~ 4 . The meaning of the symbols is explained in 
fig. | 
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A. The variable m^2 an d its generalisation 



The Cambridge m T 2 variable, proposed in [94], can be used in analyses such as [28], 
where particles are pair-produced at hadronic colliders, and decay semi-invisibly, as 
is the case in R-parity conserving SUSY. This variable is used in symmetrical, two- 
body decays of supersymmetric particles, where the LSP is unobserved, and so must 
be inferred from missing energy. 



A.l The properties of itlt2 

For events in which the decay of a heavy object produces an unseen particle, such as 



(A.l) 



one can write the Lorentz invariant 



xi 



7r rpXi 



m*_+ = m ^ + m ^ + 2 E^E* 1 cosh(Ar/) - p£ • p* : 



V 7T _Xi 



(A.2) 



x 

where p^ and p r * indicate pion and neutralino 2-vectors in the transverse plane, and 
the transverse energies are defined by 



Et=\(Pt) 2 + K and K 



Also 



is the true rapidity, so that 



rj = - log 



E + p z 
E- Pz 



(p^) 2 + m| ? 



tanhr] = p z /E , sinhr? = p z /E T , cosh i] = E / E T . 



(A.3) 



(A.4) 



(A.5) 



In a hadron collider, only the transverse components of a missing particle's mo- 
mentum can be inferred, so it is useful to define the transverse mass, 



m 2 T (p^, p* 1 ) = ml+ + m| + 2(£J£* 1 - p£ • p* 1 ) 



xi 



(A.6) 



which, because cosh(x) > 1, is less than or equal to the mass of the lightest chargino, 
with equality only when the rapidity difference between the neutralino and the pion, 
Ar)^ is zero. All other Ai] lead to mr < m ^+, so if we knew the neutralino 
momentum we could use m T to give an event by event lower bound on the lightest 
chargino mass, uit was used in this way by UA1 [95] in the measurement of the W ± 
mass. 



30 



In R-parity conserving SUSY events there are expected to be two unseen LSPs. 4 
Since only the sum of the missing transverse momentum of the two neutralinos is 
known, the variable 



m 



T2 



mm 



max{m^(p£ m^Pr ,^ 2) )} 



(A.7) 



is a lower bound on the transverse mass m-r for events where two decays of the type 
( |A.1|) occur. In (|A.7|) we have been forced to minimise over all consistent neutralino 



2-momenta. Note that is the hypothesised momentum of the ith neutralino which 
need not be equal to its true momentum. 






Figure 17: A diagram demonstrating that the minimisation over some parameter of the 
maximum of two well-behaved functions may occur either at (a) a minimum value of one 
of them, or (b) when they are equal, or (c) at the boundary of the domain. 



To find the range of values may take we first let f\ 
mUpf\ tiP) 



"4(Pr (1 \ 



4" 



We then note that the minimum over a parameter x 

= 0, 



(2) 



X 



and f 2 

of the maximum of fi(x) and f 2 {x) can occur at a local minimum, f[ 
provided fi(2)(x*) > f2(i)(%*), as shown in fig. |l7|a. Alternatively the minimum can 
occur when the functions cross one another when fx = f 2 (fig. |l7|b) or at a boundary 
(fig. |l7|c). The parameter x corresponds to the fraction of the the missing momentum 
(in one of the transverse directions) which is assigned to each half of the event. Since 
fi , f2 — > oo as x — > ±oo fig. |T7| c is not relevant to our minimisation problem. 

To see that case (a) cannot occur, consider the unconstrained minimisation over 
<fcr, of itit{Pti $t)- Using the relationship 



(A.8) 



where ffj\ = $t + m ~o, it is straightforward to show that, 

Xi 



9m T _ 9 / P7T »/•• () rr 



k = 1,2 



(A.9) 



4 Though there may also be other unseen particles - see appendix A. 2. 
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This means that at the minimum 



v? = , (A. 10) 

where we introduce the notation v-p = p^/E^, yt^ = $t/$t- where and vy 
represent the true transverse momentum and velocity of a particle, while fa and ylj 
are assigned by the minimisation. 

Using the basis (t, x, y) with the metric diag(l,-l,-l), one can write 

m 2 T = (E%\ p^) ■ (E?\ p^) , (A.ll) 

where = E£ + $t and pJp* = Pt + fa- This 1+2 dimensional Lorentz invariant 
can be evaluated in any frame boosted from the lab in the transverse plane. ( |A.10| ) 
has told us that at the unconstrained minimum the transverse velocities and 
fa are equal; a statement necessarily true in all transverse frames, including the 
special one in which both the transverse velocities (and associated momenta) are 
zero. Evaluating (|A.11|) in this frame, we find that the unconstrained minimum of 



( |A.lip then becomes (m n + m^o, 0, 0)-(m 7r +m^o, 0, 0), and we recover the expected 



result 

m ™ in = m n + m _o . (A. 12) 

We therefore conclude that the function m\, has only one stationary value and it is 
the global minimum, and is common to both sides of the event provided the same 
type of particles are emitted. Thus when fi is minimum it cannot be greater than 
/2, and so the minimisation in ( [A. 7 ) forces f\ — f%. This could of course occur when 



both fi and fi are at their global minima, in which case rriT2 takes its minimum 
value: 

m T2 min = m^ + m^o . (A. 13) 

To summarise, mx2 is the minimum of wip subject to the two constraints mJp = 
my , and pjP + pj/ = p T . The condition for the minimisation can be calculated by 
Lagrange multiplier methods, the result of which is that the velocity vectors fi T of 
the assigned neutralino momenta must satisfy 

(4 1} -v? (1) ) oc U^-vf*). (A. 14) 

To find the maximum of m,T2 over many events we note that for each event 



the minimisation will select hypothesised momenta satisfying (|A.14 ). We now note 
events can occur in which the true transverse velocities of the neutralinos were exactly 
those which were assigned by the minimisation, i. e. they can satisfy 

v| (1) =^, 4^ = ^^. (A.15) 
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These events will have both hypothesised transverse masses equal not only to each 
other but also to true transverse masses which would have been calculated if the 
neutralino momenta had been known: 

(*) rA^\ - ™« J® 



m^[pr,i>r v )=<>^r,€ } ) (A.ie) 

If events occur where, in addition to the transverse components of the neutralino 
momenta satisfying ( |A.15| ), the rapidity differences satisfy = ^t(i) anc ^ ^x?(2) = 
77^(2), then by ( |A.2| ) vtlt2 will equal the true mass of the chargino. Combining this 
with ( |A.13| ) and recalling that m<r 2 cannot be greater than the chargino mass by 
construction, we can see that the event-by event distribution of iuti can span the 
range 

m^o+m^ < m T2 < m~+ (A. 17) 

showing that rriT2 is sensitive to the m-+ — m^o = AMjj mass difference. 

The variable is equally applicable to two same-sign xt decays so rriT2 signal 
events can be defined as those having two xf Xi 7r± decays with any combination 
of charges. 

A. 2 Generalisations of 

In our analysis we also wish to make use of the leptonic decays xf Xi & v t where 
t G e, \i. We therefore generalise rriT2 to cases where more than two particles go 
undetected. 

Consider events in which a chargino is produced and then decays to Xi e v e- If 
we expand the Lorentz invariant 

(m~+) 2 = (p^o +p e +Puf (A.18) 

we obtain three mass-squared terms for each of the decay particles and three cross- 
terms. The cross-terms can each be written in the form 

Pa ■ Pb = 4 a) 4 6) cosh(Ar ?a6 ) - pP ■ pP . (A.19) 

If the neutralino and neutrino transverse momenta were individually known we could 
evaluate the transverse mass, 



m T = m~o + m e + 2 



(E e T E* - p e T ■ p*) + {E V T E* - p v T ■ p*) + {E e T E v T - p e T ■ p y T ) 

(A.20) 



where the neutrino mass is assumed to be negligible, mx will be equal to the xt 
mass in events where Ar] a b = for all pairs of e, v e , and Xv 
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Figure 18: (a) Simulations of rnxx — m x\ ^ or ^ = 2, 3, 4 using a simple phase-space 
Monte-Carlo generator program for a pair of decays q — ► Xi 9 followed by Xi ~~ > X? ^ or 
X^" — > x? e As the number of invisible particles increases the proportion of events near 
the upper limit decreases. The peak in m^z — m$o near the pion mass is explained in the 
text, (b) The distortion of ititx — m x\ w hen the LSP mass is varied by ± 10%, showing 



m- 



xi 



m^o. In this 

Xi 



that Mtx — m x\ rema i ns sensitive to the mass difference AM^ 

simulation AM^ = 0.845 GeV, m^o = 161.6 GeV, and the electron and neutrino mass 
were neglected. The normalisation is arbitrary. 



In events with two leptonic chargino decays a variable like m?^ can be defined 
as in ([A .7]) but using the three-particle definition of ttit from ( |A.20| ) and with the 
modified constraint, 



q£« + q*« + + 



(A.21) 



where the labels (1) and (2) indicate which chargino the particles were emitted from. 
We call this variable % 4 (or indeed where X is the number of undetected 

particles). 

The conditions for the minimisation required to calculate can be calculated 
just as for m^. The Euler-Lagrange equations involving 

8{mf) 2 d(m®) 2 

cLllCl 



<9q^ (i) 



show that the minimisation will select the invisible particles' momenta such that 



u 



Uy . The other E-L equations reproduce ( |A.14j) but with electrons replacing 



pions. 

This means that when calculating tuta one can replace the missing particles from 
each chargino decay with a pseudo-particle with mass equal to the sum of the masses 
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of those invisible particles and proceed as for m^. In the case of leptonic chargino 
decay the mass of the neutrino can be safely neglected in comparison to that of the 
Xi, and the constraint = Uy*' is equivalent to = (0,0). 

The generalisation to for other values of X is straightforward. The distri- 
bution of 777,7-3 is shown in fig. |T8| a for events in which one chargino decays to x\i e i v 
and another to Xv n+ ■ Unlike m T2 and m T4 it has a sharp peak at m T3 = m^o + m n . 
This occurs because the visible particles on each side of the event are different and 
so the unconstrained minimum of the values of 777,7- on each side of the event are not 
equal as they are in the case of 777,7-2 and 777,7-4: 

mm ^777^ (p^,, </£p) j = 777^ + 777^0 7^ 777 e + 777^0 = min ^m^(p^, ^P)j 

(A.22) 

Some of the events can then fall into the category shown in fig. |TT|a, producing a 
peak of events with 7777- = 777^0 + m^. 

The distribution over events of 7777-4 will have fewer entries near the upper kine- 
matic limit (7777-4 = 777-+) because when more particles go undetected an event at 
that limit must satisfy a larger number of constraints. For fully leptonic chargino 
decay, there are six constraints of the type A77 = 0, two = and finally the 



modified constraint from ( |A.14| ). This effect can be seen in fig. [L8|a for events where 



a total of two, three and four invisible particles are produced. 

A further generalisation which we do not require here might be relevant in cases 
where more than one visible particle is emitted from each mother. For such decays, 
one would sum the full 4-momenta of the visible particles from each decay as well as 
summing the masses of the invisible particles from each side of the event and proceed 
as for 7777-2. 

A. 3 Uncertainties in pV and rn^o 



The sensitivity of itltx to the estimated mass of the neutralino is shown in fig. [18|b, 
where 10% (16 GeV) errors in x result in similar fractional errors in AM^ i.e. of a 
few tens of MeV. 7777-x shows similar insensitivity to measurement uncertainties in 
the missing transverse momentum vector. This behaviour can be (at least partially) 
understood from the non-relativistic limit of 7777-2, when the proportionality in ( |A.14j ) 
becomes an equality and 



1 / \ 2 

m T2 2 - (m w + rn^o) 2 = (rnrf T - - m.^ j + O ((v T • V7-) 2 ) . 

(A.23) 



One can see that in ( |A.23 ) multiplies the missing momentum, while 777^,0 multi- 



plies the pion transverse momenta. 
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B. Simulation of particle tracks 

Part of the ATLFAST [45] software provides fast parameterised simulation of the AT- 
LAS inner detector performance. In particular it calculates the expected track re- 
construction efficiency, and smears the five track helix parameters and according to 
particle type. The parameterisation used was based on GEANT3 Monte-Carlo simula- 
tions [46] using a Kalman- filter algorithm for track reconstruction [96]. A parame- 
terisation from a large statistics Monte-Carlo sample [97], was used for the hadronic 
track smearing for this study. 

B.l Low pr electron identification 

In ATLAS the transition radiation tracker allows the identification of low-energy 
electrons because they emit more transition radiation than more massive particles 
with the same momentum. The issue of particle misidentification is not dealt with 
in the standard version of ATLFAST, and so has been implemented independently in 
our analysis. 



Vt 


0.5 


1 


2 


5 


10 


> 20 




0.01 


0.003 


0.005 


0.007 


0.008 


0.012 



Table 7: Pion misidentification efficiency at = 0.3 as a function of transverse momen- 
tum for an electron identification efficiency of 90%. Interpolations were made linearly in pr 
and logarithmically in e n . Particles with pt < 0.5 GeV are assumed not to be reconstructed 
into tracks. 



\v\ 


< 0.2 


0.3 


0.65 


0.8 


1.15 


1.4 


1.9 


2.15 


> 2.4 


p T 2 GeV 


1.2 


1.0 


0.88 


1. 


0.6 


0.4 


0.028 


0.88 


2 


p T 20 GeV 


1.1 


1.0 


1. 


2.0 


1.4 


0.78 


0.22 


1.4 


2.3 



Table 8: r\ dependent correction factor applied to e n , the hadron misidentification effi- 
ciency, for two different values of transverse momentum. Interpolations were made loga- 
rithmically in and linearly in pt and r\. Tracks beyond \r)\ = 2.5 are not reconstructed. 

Test beam performance was compared to Monte-Carlo simulations in [98] using a 
combined discriminator which considered both the number of high threshold clusters 
and the time for which the charge deposited exceeded a lower threshold. In that 
paper the pion misidentification probability was calculated as a function of transverse 
momentum at rj = 0.3 for an electron efficiency of 90% (table 0). For our simulation 
these efficiencies were extrapolated to other values of r\ by comparison with the 
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GEANT3 simulations performed in [27]. The and rj dependent correction factors 
are shown in table |8|. 



B.2 Low pt muon identification 

High energy muons are easily identified because of their high penetration. Identi- 
fication efficiency is drastically reduced when muons have an insufficient transverse 
momentum to extend their track helix into the dedicated muon detectors. The ap- 
plied identification efficiencies for muons were based on [27] and are shown in table H 



Pt 


< 2 


3 


4 


5 


6 


8 


> io 


Efficiency 





0.34 


0.66 


0.89 


0.95 


0.96 


0.98 



Table 9: Applied muon identification efficiency as a function of transverse momentum. 
Efficiencies were interpolated linearly in and pr- 
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